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Abstract 

We consider deeply inelastic scattering at very high energies in the saturation regime. The emerg- 
ing picture corresponds to the propagation of a dipole, the quark-antiquark pair, in a shock wave 
color field of the target. We use the fomalism of Wilson lines to study the evolution of dipole 
densities in energy logarithms. Our analysis results into an equation in multicolor limit which 
sums leading logs but keeps the nonlinearities up to cubic order in densities. 
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1 QCD of dense parton systems. 



A standard tool in QCD analyses of the deeply inelastic lepton-nucleon, £N — > £'X, scattering 
(DIS) for moderate values of Bjorken variable, x&, is the factorization theorems and Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equation This equation describes the 
change of the DIS observables with the change of a resolution scale, which for the process in 
question is the virtuality of the probe Q 2 . The latter serves as a microscope which allows to 
penetrate deep inside the hadron substructure and observe its constituents with transverse size 
5x± ~ Q^ 1 and longitudinal extent — . 

The DGLAP dynamics is based on the separation of the DIS amplitude into a 'hard' part 
coming from the transverse momenta k\ > fi 2 and a 'soft' part coming from low k\ < fi 2 , where 
fi 2 is a scale dividing short and long distance physics. The incoherence of these phenomena allows 
for a factorization of infrared part into a universal matrix element of a non-local composite light- 
cone operators constructed from quark and gluon fields. On the other hand, the contribution from 
hard momenta gives the coefficient functions. The factorization scale [i 2 serves as a normalization 
point for those operators. The change in /i 2 is governed by the conventional renormalization group 
equations. Taking fi 2 = Q 2 , we come to the usual result that the Q 2 dynamics of DIS cross sections 
is driven by the renormalization group equations for the light-cone operators. In terms of QCD 
perturbation theory it results from the summation of the contributions of the type {a s \nQ 2 ) n , 
etc., in momentum transfer. 

A very important property of the factorization alluded to above is that the coefficient functions 
are purely perturbative. Indeed, the effective coupling constant is determined by characteristic 
transverse momenta so that the contributions coming from large k\ > /i 2 are treatable within 
QCD perturbation theory as long as fi 2 is sufficiently large. The nonperturbative physics enters 
the game only when we lower the normalization point fi 2 down to a typical hadronic scale of order 
~ 1 GeV. The higher order terms of perturbative expansion, for both the coefficient functions 
and the anomalous dimensions of the light-cone operators, lie in the same framework of linear 
evolution and lead to corrections ~ a s , a 2 , etc. Thus, to compare experimental measurements of 
structure functions JF(xb,Q 2 ) at different Q 2 we rely only on perturbative QCD and the linear 
character of the DGLAP equations makes this comparison especially simple. 

The situation changes drastically if one is interested in the domain of small x-q- The DGLAP 
evolution leads to a strong rise of the DIS structure function 

/ i \i/2 

xbF(xbi Q 2 ) ~ exp ( const. In — lnlnQ 2 ) , (1) 

\ xb J 

at small values of Bjorken variable xb- If one bears on using the DGLAP evolution for smaller 
and smaller xb, higher loop contributions become enhanced by additional factors In ^ and the 
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perturbative expansion of the coefficient functions and anomalous dimensions breaks down calling 
for the small-xs resummation. Recall that DGLAP equation sums logs of the hard scale to all 
orders, i.e., terms of the kind (a s \nQ 2 ) n and (a s In Q 2 In . It takes only a single logarithm 
or none of the energy for a power of the coupling constant. Thus, it fails for very 1ow-Xb when 
a s In ^> 1 and these contributions have to be summed over. In perturbative QCD, the small- 
xb asymptotic behaviour is described in the leading logarithmic approximation (LLA) by the 
Balitsky-Fadin-Kuraev-Lipatov (BFKL) pomeron ||] which sums up the leading energy logarithms 

Unfortunately, the BFKL evolution, for a review, see ||, suffers from its own caveats. The 
first one is the lack of unitarity: the power behavior of the cross section due to BFKL dynamics 

x B F(x B ,Q 2 ) ~^ p , a P = l + 4iV c ^ln2, (2) 

7T 

violates the the so-called Froissart theorem stating that a cross section may grow at most as In 2 ^ 
at xb — » 0. Obviously, this result means that approximations involved in the derivation of the 
BFKL equation become inadequate and in order get the true asymptotic behaviour at small x B , 
we must go beyond the LLA. Unlike the DGLAP case, this is not a purely technical problem of 
calculating loop corrections to the kernels. There are a s corrections to the BFKL kernel but 
in addition there are unitarity corrections which go beyond the framework of the BFKL equation. 
At small a s and x-q, the latter corrections seem to dominate over the next-to- leading BFKL effects 

The second problem with the BFKL evolution is its infrared instability. We can safely apply 
perturbative QCD to the small-xs DIS if the characteristic transverse momenta of the gluons 
kj_ in the gluon ladder are large. For the first few evolution steps, one can check by an explicit 
calculation that the characteristic k\ are of the order ~ Q 2 . However, as x-q decreases, it turns 
out that the characteristic transverse momenta in the middle of the gluon ladder drift towards 
Aqcd making the application of perturbative QCD questionable. This is related to the fact that 
the operator expansion for the high-energy scattering in terms of Wilson line operators, which 
represent quarks moving with almost the velocity of light, [[| is based on the factorization in 
the rapidity, 7/ = ln^-, [0, § rather than the transverse momentum. Unlike the usual light-cone 
expansion, the high-energy expansion in Wilson operators does not admits an additional meaning 
of perturbative versus nonperturbative separation. Contrary, both the coefficient functions and 
the matrix elements have perturbative as well as nonperturbative parts. This happens because, as 
we mentioned above, the coupling constant in a scattering process is determined by the scale of the 
transverse momenta. When we use the factorization in hard (k± > /i) and soft (k± < ^1) momenta, 
we calculate the coefficient functions perturbatively, since a s (k± > fi) is small, whereas the matrix 
elements are nonperturbative. Conversely, when we factorize the amplitude in rapidity, both fast 
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and slow parts have contributions coming from regions of large and small k±. In this sense, the 
small-xe evolution in QCD is not protected from the infrared side in the same way as the DGLAP 
evolution is: in order to compare the two structure functions measured at different (small) values 
of xb the perturbative QCD may be insufficient and, in order to explain the small-xe behavior of 
structure functions, it may be necessary to take into account the interplay between the hard and 
soft pomeron. 

Both of these problems can be resolved simulteneously if, as argued in || 0, [IT], 0j the 
partons in the highly energetic nucleon reach the state of saturation: the recombination of partons 
balances the rise of the cross section due to parton emission, and the hard saturation scale Q s sets 
the scale of the effective coupling constant. Indeed, once we have a rather dense gluon system 
created by conventional parton splitting described by the linear DGLAP and BFKL evolution, the 
partons populating a given space-time volume inside the hadron start to overlap and an absorption 
competes with creation. This is expected to happen when gluons occupy the entire transverse area 



of the hadron disc S^ adr = ^-R 2 , i-e., 



Zpart (5xj_) 2 n 

where the number of partons n is proportional to the gluon density xbG(xb, Q 2 ) which overwhelms 
quarks at small xb- This is a result of parton saturation which is expected to tame the growth 
of their number. So in the case of gluon overlap, they start to interact strongly although the 
QCD coupling may well still be in the perturbative domain. This regime was addressed in the 
pioneering work by Gribov, Levin and Ryskin |§ and has resulted into a suggestion of a first 
nonlinear evolution equation which goes under their names, the GLR equation. It received an 



early discussion in Ref. [13|. A derivation of the later within double logarithmic approximation 
has been given in The question of the value of the scale at which the annihilation takes 

over the production based on the analysis and solutions || [TJ], [H| to the GLR equation has 
resulted into the aforementioned concept of saturation. 

A new era of investigations of nonlinear QCD phenomena has been been initiated by the 
experimental results from HERA on small-xe measurement of the DIS structure functions. It 
created diverse dynamical approaches to the corresponding physics which we will discuss below 
and has led a phenomenologically successful saturation models, see e.g., P"7[ . 

The most transparent picture of the underlying nonlinear effects arises in the dipole frame 



lq, 19], where the hadron carries almost the entire rapidity, and moves with almost the speed of 



light, but still the virtual photon is quite energetic. The incoming photon fluctuates into quark- 
antiquark pair, a dipole, and interacts with the target via gluon emission. Since the hadron is 
contracted due to Lorentz dilation, the dipole with accompanying radiation sees it as a color 
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Figure 1: Scattering of the photon probe off the hadron at high energies. The virtual quantum 
fluctuates into a bunch of dipoles which interact with the target. 



source of a transverse extent living on the light-cone, see Fig. [T]. This color field is created by the 
constituents of the well developed wave function of the hadron which in view of its high intensity, 
i.e., big occupation numbers, can be considered as classical. The strength of the field in the regime 
of saturation is 1/ JoT s - This can be achieved either by going to the very low x-q for a hadron 



or scattering off a nucleus [|TTJ, or both. Actually in the present paper we will consider the last 
possibility which justifies the neglect of the so-called Pomeron loops which are suppressed then by 
the atomic number of the nucleus as compared to exchange contributions. At high energies, we 
can neglect the recoil of emitted gluons so that it is legitimate to use the eikonal approximation 
20|, 2TJ|. The proper degrees of freedom for the fast particles moving along the straight trajectories 



collinear to their velocities are the infinite ordered gauge factors, mentioned earlier with respect 
to an operator approach to the high energy scattering. In this case, the color dipole is given as 
two Wilson lines stretched along the light-like direction = -^(1, 0, 1), 

U (oo, — oo; x) = Pexp (ig J dx_A + (x^, x)^j , (3) 

and separated by the transverse distance — y±^ = (0, x — y, 0), i.e., 

Af(x, y) = (0|T irtr (lA (oo, -oo; x) ifi (oo, -oo; y) - l) |0> . (4) 

It is evaluated in the external field of the hadron or nucleus alluded to above. The structure 
functions are given by a convolution of the probability for the photon to fluctuate into the quark- 
antiquark pair and the dipole cross section expressed as an integral of the dipole density over the 
impact parameter, b = | (x + y). Namely, ||22|| , 



^b,Q 2 ) ~ Jd 2 zJ dz\^ q (z,z,Q 2 )\ 2 J d 2 bM(x,y), (5) 

with z = x — y. 

In doing this procedure one introduces a separation in rapidity as discussed previously, ac- 
cording to which the gluons of high rapidity go to the impact factor, — the square of the photon 
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wave function alluded to above, — while the slow fields form the dipole. The change of this divide 
cannot change a physical observable and it is governed by an evolution equation for A/". To find 
the dipole density at higher rapidity 77 + Sr], i.e., smaller Bjorken variables, in terms of the one at 
rj we have to integrate the fast gluon modes out in the strip Srj. At low densities, this evolution 



stems from an independent branching of dipoles and results into the linear BFKL equation [23 
Saturation effects manifest themselves in breaking of this pattern when the dipoles start to 'feel' 
each other and intensively interact. Thus, there are two types of corrections to this result: First, 
radiative corrections to the photon wave functions, and second, contributions of higher, or multi- 
ple, dipole densities to the cross section. The latter arise as soon as an additional gluon is emitted 
along from one of the quarks in the pair. In the multicolor limit this quark-gluon-antiquark system 
is reduced to the pair of dipole densities both of which interact with the target field, n extra gluons 
lead to n+1 dipoles, etc., which are accompanied in the structure function T by the multi-parton 
photon wave functions, i.e., schematically, |^( ng ) 9 | 2 A/' n+1 . 

A number of studies along this line has led to a nonlinear equation which generalized GLR 
equation. It was derived first in M where an infinite set of coupled integral-differential equations 
has been given and a perturbative kernel for the first nonlinearity has been found. In Ref. |24] the 



latter was deduced from the Mueller's nonlinear equation for the generating functional of dipole 
densities |23J in the multicolor limit. In (F25] it was rederived using the direct summation of the 
fan diagrams. In |2J], [27| it was deduced from the functional equation, for the statistical weights 
of the Color Glass Condensate, a state of dense gluon matter at high energies for heavy nucleus 
p7| , 28 1 , derived in [29] and rederived in [^, [27] . The paper |3(| dealt with the infinite hierarchy 



of equations alluded to above which were summarized in a compact form of the functional Fokker- 
Plank equation and its properties has been discussed. 

The most general nonlinear evolution equation in the large- A^, limit reads 

j^±N(x, y)= J d 2 z K. x {x, z, y) {M(x, z) + Af(z, y)} + J d 2 z )C 2 (x, z, y) Af(x, z)M(z, y) 

+ J d 2 z d 2 z' K 3 (x, z, z', y) Af(x, z)M(z, z')M(z', y) + . . . . (6) 

The necessity to take the multicolor limit will be explained in the main text. Here it is sufficient 
to say that it allows the factorization of multiple Wilson line correlations into dipole densities and 
leads to a closed equation (Q). Note that the term K n starts from the a n ~ x order in perturbation 
theory. The known evolution kernels are K\ and K.2- The former, obviously, coincides with the 
BFKL evolution kernel and is given by 

Ki{x,z,y) = ^N e ( {x ~f (7) 
27T 2 I \z — x) l \z — y) 2 



-1 (5^(z -x)+ 5W(z - y)) [ dV ,„ {X *J J + O (a 2 



2V v ' y a/ H {z'-x) 2 {z'-y 



Use also J\f(x, x) = 0. The kernel for the first nonlinearity has been found in a number of studies 
to be 

^■»-») = ^ (,- ( ^ w) , + o(^)- w 

Extensive numerical and semi-analytical solutions to the generalized nonlinear equation with 



fc n >2 = performed in |3T], [25|, [3^, [33], |34], |3oj have demonstrated the desired suppression of the 
growth in the parton densities and resulted into quantitative estimates of the saturation scale 
Q s where this turnover actually takes place. In Ref. |36| the solution to the equation has been 



expressed in terms of a path integral suitable for lattice evaluations. 

Our current study is devoted to the computation of the evolution kernel /C3. Presently, we 
neglect the next-to-leading BFKL terms in JC\ and 0(a s ) corrections to the three-pomeron vertex 
/C2 and reserve their evaluation for our future study. However, one can expect them to be para- 
metrically less important that /C3 since a bulk of these 0(a s ) corrections comes from the effects 
of running coupling constant. For rather small coupling constant, the unitarity effects, which cure 
the strong rise (0), become dominant at rapidities 



2 , 1 
r] ~ In — , 

cup - 1 a s 

which are smaller than the rapidities when the running of the QCD coupling in the BFKL equation 
starts to be important. The latter was estimated to be || 

Therefore, the unitarity restoring effect dominate long before the diffusion constraints, and the 
question of unitarity restoration can be considered within the fixed coupling approximation [|37]]. 
Those readers who are not interested in technical details of the derivation, leading to the fi- 



nal result, can skip sections the most of the consequent presentation and go to section |6.7| and 
conclusions. For the rest, our paper is organized as follows. In the next section we consider the 
Regge limit of the deeply inelastic scattering amplitude and define the gluon field configuration 
of the boosted hadron which serves as a scattering source for the virtual photon. Then, in section 
|3] we compute the exact quark and gluon propagators in this shock-wave field. We use them in 
section |] for a calculation of leading and next-to-leading nonlinearities in the generalized evolution 
equation. Finally, we present a few concluding remarks. The two technical appendices give a de- 
tailed account of two-loop computation of diagrams with nonabelian vertices and two-dimensional 
Fourier transformation from the momentum to the coordinate space. 



6 



2 High-energy limit and shock- wave. 



As we already explained in detail earlier in the dipole frame the external electromagnetic probe, 
the photon, having momentum q fluctuates into a quark- ant iquark pair and passes through the 
Lorentz contracted color field of the hadron or nucleus with momentum p. Let us determine this 
field configuration in the high-energy limit. To this end we consider the path integral which defines 
the DIS amplitude 

T^(q,p) = i J d A ze^(p\T 3 ^z/2) ]u (-z/2)\p) A 

= i J d A ze tz -i(p\ J DA Dip Dip exp [i J d A y j,{z/2)j v {-z/2)\p) . (9) 

Each field configuration is weighted with the exponential of the QCD Lagrangian £ = £ cl + Cg. 
Here the classical part is 

Cd = (10) 

and the gauge fixing part will be specified later in this section complying with the form of the 
external color source which will make the treatment especially simple in the present circumstances. 
Our conventions are fixed by the following definitions for the covariant derivative = 8^ — igA^ 
and the field strength tensor G® u = d^A a v - d u A^ + gf abc ' A\A%. 

Motivated by the demand outlined in the introduction we introduce two light-like vectors 
and n* such that n 2 = n* 2 = 0, and n • n* = 1. We use the conventions for light-cone coordinates 
z + = z ■ n, z_ = z ■ n*, so that z^ = Z-n^ + z + n* + zj_^. Thus z 2 = 2z + Z- — z 2 , where z\ = —z 2 . 
The projection on the transverse plane is done with the transverse metric g^ v = g^u — n^n^ — n^n*^. 

In massless limit the target four-momentum is light-like p = n* and the virtual photon vector 
can be decomposed in Sudakov variables as follows 

% = ~x B n; - |-n,. , 

with xb = —q 2 /{2q ■ p). The small-xe limit, A = 1/xb — > oo, can be treated as a rescaling of the 
fields in the functional integral |38] 

lim T uu (q,p) = i I d 4 ze- iz — iz+ * '~ /2 

x lim (p\Tj„ (\- 1 z+/2, \z„/2, z/2) j v (-\~ 1 z + /2, -Az_/2, -z/2) \p) . (11) 

A— >0O v ' v ' 

After changing the variables in the path integral we have the amplitude 

Km T^(q,p) = i [ d 4 ze- iz -- i *+« 2 / 2 {p\Tj IM (z/2)j u (-z/2)\p) B . (12) 
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evaluated in the background gluon field 



B-(z + ,z_,z) 



B+{z + ,z_,z 



B(z + ,z_,z) 




) 



(13) 



Taking the limit A — > oo one gets from Eq. (|I3|) 



B + = 5{z.)(3{z) 



B 



B 



0. 



(14) 



Thus, the background field has the form of a shock-wave type. The field strength tensor for this 
potential 



is also localized on the light cone. 

Thus, considering the scattering of a dipole in the external shock-wave field we decompose the 
total gauge field into its background and quantum parts 



While developing the perturbation theory with respect to the quantum fields we will keep classical 
source effects exactly. To make sense out of the path integral we have to impose a gauge condition 
on the quantum field b. As it is almost obvious from these considerations, the light-like gauge 
&_ = b ■ n* = is the most appropriate choice. For this case, the addendum £ g f reads 



Since the DIS amplitude is gauge independent we take the limit £ — ► in the generating functional 
(0) which simplifies considerably the gluon propagator. 

The procedure we have just performed can be translated into the language of the Wilson 
renormalization group. We have included into the classical field the low frequency excitations of 
the gauge fields p - < ~ V x bo> while the quantum part involves fast configuration. Since 
the generating functional does not depend on the separation scale A the change in the latter is 
compensated by the renormalization group flow, i.e., in each evolution step we integrate out the 
quantum field with the momenta in the strip 




(15) 



A = B + b. 



(16) 




(17) 



Po- < p- < p - + dp 



(18) 



in order to go from the scale p _ + 5p^ to p -- 
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3 Propagators in the shock-wave background. 

As a next natural step we have to find particle propagators in the shock-wave background (|14|). 
We do this by an explicit summation of the gluon emission diagrams. For other related discussions 
of the backround gauge propagators, the reader is referred to Refs. |39], [| [50], [4_lJ, . Since we 
will use both momentum and coordinate space representations we define the Fourier transformed 
propagator according to 

G(p, p') = J d*z d 4 z' e ip - z ~ ip '- z 'G{z, z') . (19) 
And the Fourier transform of the shock-wave field is 

B;(k) = n;27r5(k^(3 a (k). (20) 

Now, we address the quark and gluon Green functions in turn. 

• Quark propagator: iS l j(x,y) = (0\Tip % (x)%l)j{i))\Q) b- 

Summing up the gluon emission diagrams with the field ([2U|) we find the following matrix repre- 
sentation for the propagator 

iS(p,p') = (2n)W 4 \p - p')iS,{p) + iS (p)A q (p,p')iS (p') , (21) 

which consists of a free term Sq and an interaction part SoA q So- Making use of the explicit form 
of the gauge potential (p0|), the amplitude reads 

oo 
N=l 



Aq(p,p') = 2vr5(p_ -p'_)7- J2 (w) 

N=l 

x / fn^l (2-)V 2 ) lp + tki- P ')jN(p,p') II P(k m ). (22) 



j=l V^" J J \ 1=1 / m=l 



Here 



* - jf / (n t) - (| - - *) (f, -) S ■ <*> 

with A = p' + — (p + J2j=x fej) /(2p-) and the standard sign-function e(x) = 6*(x) — #(— x). 
Using this result one immediately finds 

^ 9 (p,pO = 2tt5(p- -p'_) 7 _e(p_) | d 2 ze-^ p - p ') Mp_,z) - 1} , (24) 

where we introduced the notation for the color matrix depending on the external field configuration 

u i j (p-,z)=e ige( P-^ x \ (25) 
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with the gauge field being a matrix in the fundamental representation (3 l j = f3 a (t a ) 1 j. Using the 
identity 

e N (x) = 9{x) + {-l) N 9{-x) , (26) 
one reduces the propagator to the one already known ||. 

• Gluon propagator: -iGf v (x,y) = (0\Tb°(x)bl(y)\0) B . 

The free-field propagator in the light-cone gauge reads 

G °>M = p^-^ ' d ^ = ^ — - > ( 27 ) 

where the square brackets on the p_-pole stand for a particular prescription to go around it in 
perturbative computations. The main advantage of the light-like b_ = is the absence of the 
quartic BBbb interaction vertices, so that we have only triple Bbb vertex left. Moreover a simple 
analysis shows that the only relevant part of the three-gluon interaction Lagrangian is 

8C = -gf abc (d,bl) B% , (28) 

since other contributions, ~ b^B^, vanish by virtue of the orthogonality property of the light-like- 
gauge gluon propagator n*G 0iAtJ / = 0, since £> M ~ n* . 

The manipulations analogous to the one we have done previously with the fermion Green 
function give: 

(-i)G^(p,p') = (2 7 r) i 5^(p-p f )(-i)G ,^(p) + (-i)G , w (p)A 9 ,P<r(P,P')H)G ,*u(p') , (29) 
with 

A g , flu (p,p') = -2n5(p^-pi)2p_e(p_)g, u J d 2 z e"^ {p - p>) {u(p_, z) - 1} , (30) 

where we used the matrix notation for the gluon field in the adjoint representation (3 ab = if acb j3 c 
and the adjoint gauge orientation matrix is related to the ones in the fundamental representation 
of the color group via 

U ab {p„, Z) = e i9<P-)P ah {z) = 2tr f t a u(p-,z)t b u\p-, Z)\ ■ (31) 

Moreover it satisfies the following hermiticity property 

u ^\ z ) = u ba (z) , (32) 
to be used extensively later. In Eq. Q3"0| ) we have also used an obvious property 

d ^2{Pi) d w A {P2) ■ ■■d m _ im (pN-i) = <W(Pi)«W(Ptf-i) • ( 33 ) 
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Let us briefly discuss several possibilities for handling the spurious infrared l/p_-pole in the 
density matrix of the gluon propagator. Obviously, it is related to the residual gauge degree of 
freedom: one can perform a x + -independent gauge transformation which does not affect the light- 
cone gauge condition. This ambiguity can be fixed by imposing a boundary condition in x + on the 
gauge field. The vanishing of the gauge field at x + = ±00 results into the advanced or retarded 
prescription on the pole 

1 1 

(34) 



[p_] p- ± iO 

Their semi-sum results into Cauchy principal value (PV) prescription (which however does not 
follow from the path integral quantization). Next, it can be handled according to the Mandelstam- 
Leibbrandt (ML) recipe g, || 

1 M 1 4 P+ ; / 35 ) 

[p_] p- + iO ■ e(p+) P-P+ +i0' 

which puts this spurious pole on the same footing as the conventional pole in the propagator 
treated by means of causal Feynman prescription. This recipe allows for the Wick rotation in 
the Feynman integrals. The ML form of the propagator does not correspond to simple boundary 
conditions on the gauge field at x + -infinity which manifest a residual gauge freedom. However, this 
prescription has been deduced later by means of the canonical equal-time (but not the light-front) 



quantization |44| and path integral formalism by means of the Faddeev-Popov trick by changing 
gauge condition in the path integral from the temporal to the light-like gauge . The ML recipe 
results into a mild infrared behaviour as compared to strong singularities one encounters when 
the pole is hadled by means of the Cauchy principal value. 

Note, however, that in our present circumstances due to the strip restriction flI5| ) on the 
integration over the p_ component of teh momentum, the actual pole in the gluon propagator is 
not hit and all prescriptions lead to identical results. 



4 Nonlinear evolution equation. 

In order to derive the evolution equation for the dipole density, we develop a perturbation theory 
for the quantum fields integrating them out in slices of p_ momentum while keeping the external 
shock-wave field to all orders. In the Wilson line formalism, the quark traveling along the path 
x M = X-U^ + x± fl interacts with soft gluons by means of the path-ordered exponential (^|). For the 
gluon field (|i~6|) it has the form 

U (00, — 00; x) = Pexp ( ig / x) ) = U (00, 0; x) u(x)U (0, — 00; x) . (36) 
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Here on the right hand side of the equality we have used the form of the shock-wave concentrated 
on the plane x_ = so that U and u stand for the path-ordered exponentials containing the 
quantum and shock-wave fields, respectively. The quantum Wilson line are to be expanded in 
perturbation series, e.g., 

U (oo, 0; a?) = 1 + 2_,{ig) \ dx\_ I dx 2 - ■ ■ ■ / dxk~ b(xi_, x) b(x 2 -, x) . . . b(xk-, x) . 

k=1 Jo Jo Jo 

(37) 

For the hermitian conjugate we have 

, k^, \t r°° f Xl ~ r x (k-i)- 
W (oo, 0; x) = t + 2 (—ig) / dxi- dx 2 / dxk- b(xk~, x) . . . b(x 2 -, x) x) . 

JO JO JO 

(38) 

so that the unitarity property is preserved 

Z7(oo,0;aj)£/ t (oo,0;aj) = l. (39) 

Now we substitute the above expansion into the formula for the dipole 

jV(x, y) = (0\T ^tr {u (oo, -oo; x) U ] (oo, -oo; y) - t) |0) , 

and form Wick contractions. Note, that without adhering to the large iV c approximation our 
evolution equation will not be closed but rather it will involve the higher correlations of Wilson 
lines, e.g., 

jV {2) (x,z,y) = (0|r— tr(W(oo,-oo;aj)Z/t(oo,-oo;«)-l) 

x irtr (U (oo, -oo; z) (oo, -oo; y) - t) |0> , (40) 

etc. However, the multicolor limit will give a possibility to reduce all higher order correlation to 
a simple product of dipole densities, see e.g., Eq. fl5l] ) below. 



5 Leading nonlinearities. 

To start with let us recapitulate the computation of the BFKL part and leading nonlinearities of 
the generalized equation (§). To this end let us compute the diagrams given in Fig. 0. For these 
contributions we only need the ++-projection of the shock-wave propagator which reads 

p z + tU \p-\ 

-2n5(v -v') - PP'f d 2 ze~ iz < p ~ pf) 

° [P - P -\p2 + i0)(p>2 + i0) [p_] J 6 

x {%_) (u(z) - 1) - 0(-p.) (u\z) - l)} . (41) 
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We can cast the free propagator to the form when the particle propagate to an intermediate 
point z before it reaches its final destination, a form we have for the interaction part with the 
shock-wave background. To achieve this, we integrate, in the Fourier transform of the free light- 
cone-gauge propagator, over the p + -component, then insert the unity 1 = J d 2 p'5^ 2 \p — p'), and 
use the integral representation of the ^-function. These manipulations lead to the result 

W-^-Zl^e-^-^-^ ,42) 
/ <// '~ 1 {6(x- -y_)9(p-) -%_ -x_)6(-p^)}e- ip - {x - y) + 

d2 P Jp-( X -z)- lp y(2p.)x^ f d2 P' -ip'iy-z)+ip' 2 /(2p-)y- P P' 



2 J 2tt \p_ 



J (2vr) 2 J (2vr) 2 (p_) 2 

-6(x- - y-)S®(x -y)J ^_i_ e -^<*-*>+ » ( 43 ) 

where we have used a p + -independent regularization of the spurious l/[p_]-pole, e.g., the ad- 
vanced/retarded/PV prescription, not the ML which does depend on p + . However, since p- 
integration is resctricted this limitation is irrelevant. In Eq. (42) the last term originates from the 



'Coulomb force' 46 



Note that if the integration over p_ would be unrestricted then one will use a property of the 
propagator with ML prescription which follows from its causality and conclude that 

G , ++ (x + = 0,x,x.) = L 0, (44) 

which holds when one performs the p- integral in the complex plane and notices that due to the ML 
prescription both, Feynman and spurious, poles lie on the same side of the real p_ axis. This would 
lead in turn to the inability to reproduce the BFKL equation. This will hold even if one would use 
the regularization by means of displacing the path from the light-like path adding a deviation 5x + 
since the corresponding contributions are finite, they do not contain double logarithmic divergence 
as in the principal value prescription. 

Note that since we have the integration over p_ in the restricted domain flT8f ), 

rpo-+&p- dp_ ^ fp -+5p_\ ^ 1 
Jpo- P- \ Po- J %b ' 

we never hit the 'spurious' pole in the gluon propagator in the one-loop diagrams since the 'ghost' 
part of the propagator is concentrated at p_ =0, namely, 

P+ -PV— - i7ce(p + )5{p^) . 



P-P+ + iO p- 

Therefore, to the leading order in the coupling constant all prescriptions on the p_-pole are 
equivalent, as we already mentioned above. 
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(a) (6) 



Figure 2: Leading order diagrams which result into BFKL equation and the first nonlinear 
correction to the generalized nonlinear equation. 

Adding the free and interaction pieces together we get 

G ++ (x,y) = l||^J-{0(x_ -y_)0(p_) -x_)e(-p_)}e-*-t-»)+ (45) 

X f d 2 Z f c iP-(*-z)-iP 2 /(2P-)*- f c -ip'-(v-z)+tP ,2 /(2p-)y- Pj 



(2tt) 2 7 (2tt) 2 (p_) 2 

x + 0(-a;_)0(-y_) + 0(z_)0(-y_)ii(;s) + 9(-x_)6(y_)u\z)} 

-8{x_ - y_) [ ^ JL e -*-t"-*)+ / ^JL^^v) l — . (46) 
V y J J 2n [p_] 7 (2vr) 2 p_ v ; 

Since the computation of the one-loop graphs is trivial we just mention that we have used in 
the derivation the Euclidean <i-dimensional transverse space Fourier transformation, 

rfrfp e ip* i r(d/2-m) 1 

7 (2vr) d p 2m ~ 2 2m n d / 2 r(m) z d ~ 2m ' ^ 

and the following identity 

2u ab (z)ti [t a u(x)t h u\y)] = ti[u(z)u\y))ti[u(x)u\z)] - — tr [u{x)u\y]} , (48) 

stemming from the color Fiertz transformation. 

The contribution of the exchange-type diagrams (with one of them sampled in Fig. ^ (a)) reads 

. r/ \ a s , T , 1 /" ,9 / „\ (z _ a?) ■ ( z _ v) 
A^,y) = ^iV c ln-/dM-2) ( l z _^_^ 2 

x |jV(a;, *) + AT(z, y) - JV(aj, y) + JV (2) (aj, z, y)} . (49) 

The contribution of the diagram |] (6) (and analogous one with the self-energy attached to the 
other eikonal line) is 

ex. . 1 f ,o f 1 1 1 



x B J {{z-x)^ \z-yy) 

x,z) +M(z,y) -AT(x,y) + Af {2) (x, z,y)} . (50) 
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The 'Coulomb' piece does not produce the logarithmic integral at leading order and, moreover, 
can be completely eliminated at this order by a regularization of the Wilson line enforcing the 
ordering of emitted gluons by a infmitezimal cutoff, so that e.g., £2 < %i — 0+ in Eq. (|38|) . 

In these computations no multicolor approximation has been involved. However, to produce 
the closed equation for dipole densities we have to adhere to the large iV c limit, 

JV (2) (x, z,y)= AT(x, z)M{z, y) . (51) 



Summing up the expressions (|49| ) and ( p0|) together with the approximation fl5"ID we reduce 
the result to the nonlinear equation 



x {Af(x, z) + Af(z, y) - M(x, y) + J\f(x, z)M(z, y)} , (52) 
which is obviously Eq. (^) with /Ci j2 given by Eqs. (0) and @, respectively. 

6 Next-to-leading nonlinearities. 

Now we are in a position to address the next-to-leading nonlinearities to the generalized nonlinear 
equation. To this end we have to evaluate the graphs shown in Fig. ^. In their computation we 
will keep only the contributions which are bilinear in u ab stemming from the gluon propagators in 
the external field and omit the color suppressed terms after application of Fiertz identities. Thus 
we will not endeavor the calculation of the radiative corrections to the leading results which arise 
from the same diagrams but with only one or none internal gluons intersecting the shock-wave. 

6.1 Preliminaries. 

Beyond leading order all components off the gluon propagator are relevant. Therefore, we find it 
instructive to use the Sudakov decomposition for momentum and decompose the density matrix 
into the part orthogonal to the vector n* and the rest, 



dpAP) = d^ip) - 2—n^n u = - £ - 2—n^n u . (53) 



Obviously, d^ip, p + — 0) = d^ v (p). Note that the last term is relevant for the free propagator 
only since it vanishes in the interacting piece as the d's are contracted in an index d fip (p)d pu (p') = 

Similarly to the treatment of ++ component of the gauge propagator in the preceding section, 
we find that the most convenient form of the latter is to integrate out the '+'-mornenta in the 



15 








(i) 



(0 



Figure 3: Samples of diagrams of different topologies inducing the next-to-leading nonlinearity, 
cubic in dipole densities, in the generalized nonlinear equation. 

d^-part. This results into the on-shell condition for the virtuality, p + = p 2 /(2p_). In this way, we 
have for the free propagator, 

GV(* -y) = ~l d p- {0(x- - y-)0(p-) - %- - *_)0(-p_)} e~^M+ (54) 



(2tt) 2 p_ m ^ " ^ (2vr) 4 p 2 + z0[p_] 

and for the part interacting with the shock-wave background 

Gi u (x,y) = - l -j^e-^-y^jd 2 z{e{x„)6{-y^)6{ P ^u{z)-6 

X / ^e^-»)-^/«*->- / ^e-^-^ z )^V(2 P -),-l d ± (p) ^ (p / ) . (55) 

Using these rules together with the conventional expressions for the QCD interaction vertices 
we are in a position to compute the graphs displayed in Fig. fj. We divide the diagrams into 
graphs of different topologies, which we name self-energy (a-f), vertex (g-i) and box (j-l) types. 
Although this division is ambiguous we stick to it for the rest of this section. 
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6.2 Self-energy topology. 

The calculation of the gluon bubble in the external field leads to the most involved algebra. 
Leaving technical details of this exercise to Appendix we give here just the final result for the 
diagram [5] (a) which reads 



J\T(x,y) = -2^ £du Jy- J ' d 2 z d 2 z' ti {{t a u(x)t d + t d u(x)t a ) u\y) } f abc f 
d 2 p 1 d 2 p[ d 2 p 2 d 2 p' 2 Q i {Pi+P2)- x - i (p'i+p'2)-y- i (Pi-p'& z - it ~P2-p' 2 )- z ' 



def u be (z)u cf (z') 



(2VT)2 (27r)2 (2VT)2 (27r)2 (pi + p2 f (p / + p > 2 f {up 2 + ^ ^2 + ^, 

x L(l +u){ (pi - pi) p' 2 ■ p' + (p[ 2 - p' 2 2 ) p 2 ■ p + 2p 2 ■ p' p' 2 ■ p} 

+u(l + u){ {pi - p\) p\ ■ p' + (p' 2 - p' 2 ) P 1 p + 2p 1 p> p[- p} 
- 2(1 + uu)[p 2 p' p[ p + p 1 p' p' 2 ■ p} - 2(1 - uu) p ■ p' (pj ■ p 2 + p 2 ■ pi) 

+ —p p' {u 2 (l + u 2 )p 2 ■ p' 2 + u 2 {\ + u 2 ) Pl ■ p[\ \ . 



(56) 



Here 



P = Pl + f>2 , P' = P'l + P 2 • ( 57 ) 

The diagram in Fig. ^| (b) is easily obtained from (a) by keeping only the first term in the color 
trace, changing the sign of the whole contribution, and by identifying y — > x in the exponential 
in the integrand, i.e., 

_ e i{Pl+P2)- x - i {p' 1 +P'2)- x -' i {Pl-Pl)- z - i {P2-p'2)- z ' _ (5g) 

When the gluon self-energy is attached to another Wilson line we keep the second contribution in 
the color trace, set x — > y, and multiply the result by the minus sign, i.e., 

_ e *(Pi+P2)-y-*(pi+P2)-y- i (Pi-pi)- z - i (P2-P2)- z ' . (59) 

The two-loop quark self-energy diagrams |3] (c), (<f), (e) and (/) are computed along the same 
line. The diagram (c) gives 

N{x,y) = 16^ C du (— I d 2 z d 2 z' tr {t a t b u(x)t c t d u\y)\ u ad (z)u bc (z') (60) 
N c Jo J p- J L > 

I d 2 p l d 2 p\ d 2 p 2 d 2 p' 2 Pl ■ p\ p 2 ■ p' 2 

J (2tt) 2 (2tt) 2 (2tt)' (2tt)2 {U p 2 + up 2 ) {up' 2 + up' 2 ) 

x i e i(p 1 -p' i y( X -z)+i(p 2 -pi 2 y( X -z')'^ 1 + e i(p 1 -p' i y(y-z)+i(p 2 -p' 2 )iy-z')'^ 1 1 

I «p? Pl 2 u pl p'2 2 J ' 

The diagram (d) can be deduced from (g), to be computed below, and reads 

JV(x,y) = 8^ I 1 du (— ( d 2 z d 2 z'u ad (z)u be (z') (61) 
N r Jo J V- J 



17 



d 2 Pl d 2 p[ d 2 p 2 d 2 p' 2 e -i( Pl -p' 1 yz-i(p 2 -p^-z' 



(2tt) 2 (27T)2 (27T)2 ( 2 7t) 2 ( p / + p / )2 ( - p 2 + up 2) ^2 + ^ 

x [uup 1 ■ p 2 (p[ 2 - p' 2 ) -2up x - p\ p 2 ■ (p'i + p' 2 ) +2up 2 - p' 2 p l ■ (pi + p' 2 )} 
lif cde tr{t a t b u{x)t c u\y)} |J_ e i (Pi+P2-pi-p^)^ _ J^giCPx+Pa-pi-pi)^ 



+ if cab ti\t c u(x)t d t e uHy)\ ( J—^Pi+Pi-P'.-P'J-y _ JL 
L J \up{ upi 



^(Pi+P-z-p'i-p'zYk 



Finally, the diagrams § (e) result into 
v 2 r x r 



N{x,y) = 16^ ( du ( ^- j ' d 2 z d 2 z' tr [t a t b u{x)t c t d u ] {y)} u ac {z)u bd {z') (62) 

J.\q '-'0 J *p — J 

d 2 Pi d 2 p[ d 2 p 2 d 2 p' 2 Pl ■ p[ p 2 ■ p' 2 



x 



(2vr) 2 (2vr) 2 (2ir) 2 (2tt) 2 + upl) ^ p ,2 + up ,2 

i e *(Pl-Pi)-( a; - z )+«(P2-P2)-( 3; " 2: ') - _)_ e i (Pl~P'l)-(y- 2: )+ i (P2-P2)-(j'- 2: ') _ 

I Pl Pl Pl 



, 2 [ ' 
Pl J 



and (/) is 

„2 f i f 



Af(x, y) = 16^- / du f^L I ' d 2 z d 2 z' tr {t a u(x)t h t c u\y)t d ) u ab (z)u dc (z') (63) 

/• rf 2 Pl d 2 pj c? 2 p 2 dy 2 i^-ptvfa-^+i^-rtvfi,-^ 1 Pi • gj P2 ■ p'2 

J (2tt) 2 (2tt) 2 (2tt) 2 (2tt) 2 uu p 2 p 2 p' 2 p' 2 ' 

The color algebra will be done in section p.5|. 



6.3 Vertex topology. 

The diagrams having vertex topology are dealt with the same technique. We refer the reader the 
Appendix [X] for the computation of the diagrams ^| (g) and present only the final result here. 
Namely, we have 



Mx,y) = -8^- C du [ d 2 z d 2 z'u ad (z)u be (z') (64) 

N c Jo J p- J 

r d 2 Pl d 2 p[ d 2 p 2 d 2 p' 2 e -i( Pl - P ' 1 yz-t(p 2 -p' 2 yz' 
X J (2vr) 2 (2vr) 2 (2vr) 2 (2vr) 2 fa + p ' 2 f ( fip 2 + Mp 2 } ^2 + ^ 

x [uup 1 ■ p 2 (pf - p' 2 2 ) - 2u Pl ■ pi p 2 • (pi + p' 2 ) + 2m p 2 ■ p' 2 p x ■ (pi + p' 2 ) 
x ^/^^{^^(aj)^^)} (^'^^^'^^^ " ^ e i(p 1 +P 2 y y -i(p' 1 + P ' 2 y^j 

+ if cab tr \t c u{x)t d t e uHy)\ \ J^^Pi+Pzyy-^+p^ _ J_ e i(pi+P2)^-i(p'i+p^-y ) 
1 > \ up\ up 2 I 
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for 0(0), 

„2 r l ,. ^ 



N{x,y) = 8^ / du [ ^ [ d 2 z d 2 z'u ad {z)u be {z') (65) 

/ 6o ^r{it(a;)t d * e ut(t/)t c } +i/ dce tr{ M (x)t c n t (y)t 6 t a } 1 



x 



d 2 Pi ^Vi d 2 p 2 d 2 p' 2 e -i(P 1 -P , 1 )^-i(P 2 -p^)^' 



(2vr) 2 (2tt) 2 (2tt) 2 ( 2 7r) 2 p?p 2 ^ + ^ + ^ a 

x ( e i( - Pl ~ p ' 1 ~ P2 ' ) ' 3:+ip2y — e ^ P2_p 'i _p 2)' :y+ip i a; ^ 
x — {wupi ■ p 2 (pi 2 - p' 2 2 ) - 2«p 1 ■ pi p 2 ■ (pi + p' 2 ) + 2up 2 • p' 2 p x • (pi + p' 2 )} , 

for FJ (h), and 

7V(ac, y) = -16^ J* du j ' ^- J d 2 z d 2 z' tr {t a t b u{x)f t d u\y)} u bd (z)u ac (z') (66) 

f^fp^fp^fp^ nr^yz-np^yz' 
1 (2vr) 2 (2tt) 2 (2tt) 2 (2tt) 2 



x 



eiCPi+Pa-p^ x-ipi-j/l Pi ■ Pi P2 ■ P2 

{ u p\ p'i P2 2 (up 2 + upl) 

+ e i(p 1 -p' 1 -p' 2 ) a ! +ip 2 yl Pi ■ Pi P2 ■ P2 

u pj p\ pi 2 (wpi 2 + up' 2 



_j_ e i Pr c C +i(p 2 -p' 1 -p' 2 )-3/I Pi ■ Pi P 2 ■ P 2 

w pi p 2 p 2 2 (V1 2 + rf) 

_j_ e -iP2-a=+i(Pi+P2-Pi)-X/I Pi • Pi P2 • P2 \ 

u p\ pi 2 p' 2 2 (up\ + wpl) J ' 

for FJ (i), respectively. Note that the last contribution is nonleading in multicolor limit according 
to Eq. (|?2]). 

6.4 Box topology. 

Finally, the calculation of diagrams having the box topology is straightforward and we get the 
result for FJ (j) 

v 2 r l r c [p_ 



AT(x,y) = 16^ f du f — j d 2 z d 2 z' tr {t a u(x)t b t c u\y)t d ) u ac {z)u db {z') (67) 

f d 2 Pi d 2 p'i d 2 p 2 (Pp^ i(p 1+P2 ). a! -i(p' 1 +p 2 ).y-i(p 1 -p' 1 ). z -i(p 2 -p 2 ). z ' 1 P1P1P2 P2 
J (2vr) 2 (2tt) 2 (2tt) 2 (2tt) 2 uu p\ p\ p' 2 p' 2 ' 



Here we have used the property fl52"|). As seen from the topology of the graph, this diagram is 
non-planar as a result it is supressed in N c as compared to leading ones, see Eq. (|72|) . For the 
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diagrams of the type in Fig. ^] (k) we have 

M{x,y) = 16^ f du f d 2 z d 2 z' ti [t a t b u(x)t c t d u\y)) 

d 2 p 1 d 2 p\ d 2 p 2 d 2 p' 2 e i ( p i +p 2)' a; - i (p'i+P2)-y~ i (Pi~p'i)- 2: ~ i (P2-P2)- z ' 



x 



/ 



(2tt) 2 (2tt) 2 (2tt) 2 (2tt) 2 ( Mp 2 + Mp 2 } ^ p /2 + Mp ,2 



x Pl ■ P P 2 • p'A u a \z)u bd {z') + + ^(*KV)^2 • (68) 

I V m pipi m p1p 2 / pip 2 J 

The color algebra of the first term in the curly brackets is handled according to the Eq. (|72|) and 
is suppressed, while the second one is reduced via Eq. ( |73|) . Last but not least, the diagrams ^| (/) 
lead to 

„2 r l r dp 

V- 

d 2 Pl d 2 p[ d 2 p 2 d 2 p' 2 ^^y^-^y 
(2tt) 2 (2tt) 2 (2tt) 2 (2tt) 2 



A/"(aj,i/) = -32^ f du f— j d 2 z d 2 z' tr {t a t h u{x)t c t d u\y)\ u ad (z)u bc (z') (69) 



\ upj p' 2 p' 2 (up 2 + upl) 



4. e i(Pi-Pi-P , 2 )-V+*P2^I Pl ■ Pl P2 • P2 

u Pi Pl P2 (up' 2 + 1 

Having computed all contributions we are ready to discuss their color properties. 



6.5 Color algebra. 

To figure out which contributions are actually suppressed in the large N c limit we simply use an 
obvious equation t a = -i^-J abc t b t c and the color Fiertz identity (t a ) i j (t a ) k t = \5 l i5 k j - 2 ^S i j S k i 
a number of times until we get rid of all Gell-Mann matrices in between the it-matrices. This 
procedure leads to the following reduction formula for the structure of diagrams (a) and (b) 

4tr {t a u{x)t d u\y)\ f abc f def u be {z)u cf {z') = tr [u(x)u\z)) tr [u(z)u\z')} tr [u(z')u\y)} 

+ tr {u(x)u\z')} tr {u{z')u\z)} tr {u(z)u\y)} 
-tr {u(x)u ] (z)u(z')u ] (y)u(z)u ] (z')} - tr [u(x)u ] (z')u(z)u ] (y)u(z')u ] (z)} . (70) 

Obviously, the last two terms in this equality are 1/iV 2 power suppressed as compared to the first 
two. 

The reduction of other diagrams is accomplished by means of the results, 

Aif cab ti {t c u{x)t d t e u\y)}u ad {z)u be {z') = tr {u{x)u\z)u{z')u\y)u{z)u\z')) (71) 

— tr $yu(x)u\z)^ tr {u(z)v)(z')} tr {u{z')u\y)\ , 
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for diagrams (d), (g) and (h), 



4tr {t a u(x)t b t c u\y)t d } u ac (z)u db (z') = tr {u(x)u\z')u(z)u\y)u(z')u\z)} (72) 
~Y tT {u(x)v!(z)}tr{u(z)v!{y)} - — tr [u(x)v! (z')} tr [u{z')v! (y)} + — tr {^^(j/)} , 

for (e), (z), (J) and (A;), and finally 

4tr {t a t b u{x)t c t d u\y)) u ad (z)u bc (z') = tr {u(x)u\z'))tr {u(z')u\z))tr {u(z)u\y)) (73) 
—ffto {u{x)u\z)} tr [u(z)u\y)} - —tr{u(x)u^(z')}tr{u(z')u^(y)} + —tr{u(x)u'<(y)} , 

for (c), (/), (fc) and (I). 

In our consequent discussion we limit ourselves to the consideration of the multicolor limit. 
Therefore, we suppress all 1/N C effects in the above equations, i.e., we keep the first two terms in 
Eq. (|70|) , the second one in Eq. fl7ll) , and the first in (|73|). 

6.6 Subtraction of multi-Regge kinematics. 

An immediate feature which is transparent in the expressions we have derived for contributions 
of particular graphs is the presence of the double-logarithmic situation, (/ dp-/p-) (J du/u). Ob- 
viously, this is an expected result and it corresponds to the iteration of the non-linear term of the 
leading order equation K 2 N ® JC 2 N 2 + /C 2 A/" 2 <g> K 2 N: 

x B ) J d Zd Z (x - z) 2 (z - z') 2 (z' - y) 2 ' 

(74) 

It has to be subtracted. The present case is similar to the one encountered in the computation of 
next-to-leading logarithmic corrections to the BFKL kernel [17 and was handled by cutting off 
the ^-integral 

pi— 6 

du — > / du , (75) 



m*. v) = 2 (^Ailn -I)' / f, ^ ^, z>W(z>,y) . 



s 



which corresponds to a cutoff in the invariant mass of the produced parton system, and omitting 
the term coeff. ■ In 1/8. 

The double logarithmic part of the diagrams we have just computed is 



A/>, y) = -2 Q^#c) 2 In In ~ J d 2 z d 2 z' (t J {a ^ Af(x, z)N{z, z')M{z\ y) 
with given by the expressions on the diagram-by-diagram basis 



(76) 



7 1 { (x - z') ■ (y - z>) (x-z)-{y-z) \ 

{a) {z-z'f \(x-z') 2 (y-z') 2 (x - z) 2 (y - z) 2 j ' 1 ' 
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J(b) 

J(c) 

J(d) 



1111 

+ ~, 777T + ~, 77T + 



2 (z' - z) 2 \ (x - z) 2 (x-z') 2 (y-z) 2 (y - z') 2 
11 11 



2 (x-z) 2 (x-z 1 ) 2 2(y - z) 2 (y - z>) 2 



1 (z> - y) ■ (z' - z) + I 



\z — X • (Z — z 



(y - z) 2 (z' - y) 2 (z' - z) 2 (z' - x) 2 (z - x) 2 (z - z') 2 
1 



(x — z) 2 (y — z') 2 ' 

(x-z) - [z! - z)(x- z') ■ (y - z') (x - z) ■ (y - z) (z - z') ■ (y - zf) 



x - z) 2 (z' - z) 2 (x - z') 2 (y - z') 2 (x - z) 2 (y - z) 2 (z - z') 2 (y - z') 2 
1 (z' - y) ■ (z' - z) 1 (z-x)-(z- z') 



(x — z) 2 (z' — y) 2 (z' — z) 2 (z' — y) 2 (z — x) 2 (z — z') 2 ' 



J(i) 



I (z 1 - x) ■ (z' -y) 1 (z-x) ■ (z- y) 



(x — z) 2 (z' — x) 2 (z' — y) 2 



z ' ~ y) 2 ( z — x) 2 (z — y) 2 



In order to compute these Fourier transforms it was enough to use Eq. ([47]). Summing up these 
contributions we obtain the required result 

(x - y) 2 



6.7 Evolution kernel /G. 



(x - z) 2 (z - z') 2 (z' - yf 



(7S 



After the multi-Regge kinematics being subtracted, the remainder defines the evolution kernel in 
the generalized nonlinear equation of the term trilinear in dipole densities, i.e. /C3. 

The afore mentioned subtraction of multi-Regge kinematics corresponds to the regularization 
of the singularities in u a la +- prescription. Therefore, in order to extract the quasi-multi-Regge 
region one has to make the substitution 



1 

u 



1 ^ dv 1 

--5(u) / — , - 
u Jo V u 



1 f 1 dv 

--5(u) / — , 
u Jo V 



(79) 



in Eqs. 

We transform the results into the coordinate representation, see Appendix [B] for details. Ex- 
tracting the factor of coupling constant from the kernel 



K, 3 (x,z,z',y) 



2T 2 ' 



N c {K 3 (x : z, z', y) + K 3 (y, z\ z,x)} , 



(80) 



with 



K 3 = J2 K : 



3(a) 
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and where we have on the diagram-by-diagram basis for large- N c part, 

K 3 (a) 



2(z - z'Y 

+ {2(x-z)-(y- z) y - ; \ 1 + 2{x - z') ■ (y - z'Y 



[x — z) 2 " (x — z'Y 



+ ((x -z).(x- z>) ((x - zf + (x- z') 2 ) -(x- z)\x - z') 2 ) [I _ ± JV_ g 

- (x- z'Y(y - z) ■ (z' -z)-(x- zY(y - z') ■ (z - z 1 ) + 2(z - z'Y 

- 4(a; - z) ■(y-z){x-z)-{y- z') + 4(aj - «') • (y - z') (a; - z) • (y - z) 

- 2(z - z'Y ((x - Z ).(y-z) + (x- z') ■ (y - z')) j 

£(x) 



/\2 



X 



(z - z') 4 ((x - zY(y - z'Y ~(x- z') 2 (y - z) 2 ) ' 
1 



= - ( x-zY( Z - z r( Z >-*A {z ~ x)iz> - x)+2 {{z ~ x) ■ {z ' ~ x)) 

+ ((z - a;) 2 + (z' - xY) ((z - z') 2 - (z - x) • (z' - x)) | 

+ g(z - z'Y + (*-*)■ (z> -*)- {x }* x -g*_ x)2 ((« - xf + V - X)' 

x 



#3(c) = 
#3(d) = 



z-z') 4 ^-^) 2 -^' -^O 2 ) ' 
1 + £(a;) 



:ilr ' (z - xY(z' - x) 2 ' 

(x-z)-(x- z') + 2 (z - z') 2 



+ 2 



(a;-z) 2 (z-z') 2 (:z-z') 2 

(a; — z) • (z' — z) (x — z) ■ (a; — z') ) £(a?) 



(aj-z') 2 (a;-z) 2 -( : z;-z / ) 2 j(a;-z) 2 (z-z / ) 2 ' 

\, ^ , txi x - z'Y , , .Ay — zY — (y - z') 

K m = -Uy-z)-{y- z'Y- + (x - z • (x - z')f i—^- i 

vy; \ (y — z'Y yx — zY — yx — z') 



+ 2 [x-zYiy-z'Y ( {V ~ Z ' nZ ~ Z ' ] ' ( " " Z ' ] + (X - Zf{Z ~ Z) - {y ~ Z 

+ 2 t-%{y-z'Y ( (y - z ' )2{z ' -*)■(«-*) + (*- *0 V -*)-(y 

£(x) 

x 



(z - z'Y((x - zY(y - z'Y -{x- z'Y(y - z) 2 ) ' 

x — z') ■ (y — z')(x — z) ■ (z' — z) 



K m = - \(x - z') ■ (y - z') + (z - z') ■ (y - z') + 2 



X 



(x - z' f 
£{x) 



(z - z'Y(x - zY(y - z'Y ' 
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K = £ [z -x)-{z - y) [z-x)-{z- y) 

3(fc) (z - x) 2 (z' - yY {x - z)\y - z'Y - (x - z'Y{y - zf K ' 



2 (z' - x) ■ (z' - y 

m ~ Jx~ 

with 



K W) = 7 wr-, ^K x ) 

(x — z) 2 [z — x) 2 (z' — y) 2 



(x) = In ■ 



[z - x) 2 



z' — x) 2 



The sum of these contributions results into the total kernel /C3. The log- free piece of which 
simplifies considerably and reads — | (D^AO ( z — z')~ 4 . On the other we did not observe an 
essential simplification for the logarithmic part. At this point it is timely to say that the total 
kernel satisfies an important property. It vanishes in the limit y — » x, 

ICz(x, z, z\ x) = , 



reflecting the unitarity property of the Wilson lines (|39|). Namely, as y — » x: the diagrams (a) 
cancels with (6), and (d) with (g), respectively. Next, the log-free term of (c) sums to zero with 
(k). Finally, the remaining contributions vanish independently. 



7 Conclusions. 

In the present paper we have developed a formalism which allows to evaluate successively the 
nonlinearities in the generalized evolution equation for the dipole densities. As a demonstration 
of our machinery we have calculated the kernel /C3 which enters with the cubic nonlinearity in the 
above equation. Presently, we have not discussed the question of inclusion of the running of the 
coupling constant into our formalism since it runs beyond the scope of this paper and requires a 
computation of radiative corrections. This will be done elsewhere. 

An obvious continuation of our analysis is to perform a (numerical) study of the evolution 
equation keeping the /C3 contribution and observe how this affects saturation phenomena. 

Other major problems for further research include: (i) A derivation of the above equation from 
the Mueller's dipole model by computing the radiative corrections to the dipole decay kernel, (ii) 
A computation oil/N c corrections to our result, (iii) A study of the effects due to derealization 
of the color source from the light cone. Recall that the latter has the shock-wave form flUp only 
in the asymptotic limit — ► 0. 
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A Calculation of diagrams with nonabelian vertices. 

Let us note first that the three-gluon vertices will not contain the '+ '-components of momenta 
after contraction with gluon propagators since the latter are orthogonal to n*, i.e., 

-?-r = o 

Since the three-gluon vertices lie on different sides with respect to the shock wave, we have z_ > 0, 
z'_ < 0, for x_ > and y_ < 0. Then, as can be seen from the explicit form of Eq. ([55]) the ' — '- 
component of the shock- wave propagators traveling through the external field is positive pj_ > 0, 
% — 1,2. As a result of momentum conservation for the '—'-components in the vertices, the 
momentum of the free external propagator = pi_ + p2~ > 0. Therefore, one finds that the 
n*n*-part of the free propagator, see Eq. (|54]) , does not contribute since the poles in p 3+ lie on 
the same side of the imaginary axis. Namely, 



r dx. r dz_ z 00 rfp3+ ^_ e -^3 + ^-+ J ( P3+ - P ?/(2 P1 -)- P i/(2 P2 _)),„ (A1) 

Jo Jo J-oo p\ + i0 

oo ] 

dP3+l TJTT; \ TTJT, \ . *7v\~~7 T—TTT = ■ 



2p 3 - -/-oo (p 3+ - p?/(2 Pl _) - pi/(2p 2 _) + i0) (p 3+ - P§/(2p 3 -) + iO) 
Same result holds for another free propagator connnecting the loop to the Wilson line. 

Thus, the only nonvanishing contribution is generated by the d x part of the free propagators 
and we have finally 

J\f(x,y) = -2^- £duuu J dp^p^J d 2 z d 2 z' tr {t a u{x)t d u\y)} f abc f def u be (z)u cf (z') 

f d 2 p 1 d 2 p\ d 2 p 2 d 2 p' 2 i( Pl+P2 ). a! -i(p' 1+? / 2 ). 1/ -i(p 1 -p' 1 ). 2; -i(p 2 -p' 2 ). z ' 
J (2tt) 2 (2vr) 2 (2vr) 2 (2vr) 2 

v Dl(PuP2,p' 1 ,P 2 ) ft r>\ 

(Pi + P2) (Pi + P 2 ) (uPi + up 2 2 ) [up\ + up' 2 J 

where we have used the substitution pi_ = wp_ and pa_ = (1 — w)p_ in order to reduce two 
integrals w.r.t. pj_ by 

/•oo roo rco rl 

/ dpi_ \ dpi- = / dp-p- / dit . (A. 3) 

jo jo jo jo 



Here 



Di(pi,p 2 ,p 1 ,p 2 ) = r^fa, -p! - p 2 ,P2)r pxa (p' 1 , -p[ -p' 2 ,p' 2 ) 

x d^ (p! + pa) (pO djL, (pi) <# fl (pa) 4 CT (p 2 ) d^ + (p[ + p' 2 ) , (A.4) 
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where we have introduced the notation for the three-gluon vertex 



r^p(Pl,P2,Ps) = (Pi - P2)p9nv + (P2 - Pz)n9vp + (Pa - Px) v g P n ■ (A.5) 

A simple computation leads to the result 

p 2 _ (uu) D 1 = u(l + u){ (pi - pfj p' 2 ■ p + (p[ 2 - p' 2 2 ) p 2 • p+ 2p 2 • p' p' 2 ■ p} 
+ «(1 + tZ){ (p\ - pi) p[ ■ p' + (p' 2 2 - p\ 2 ) Pl p + 2p 1 p' p[- p} 
- 2(1 + uu){p 2 ■ p' p\ ■ p + p l ■ p p 2 ■ p} - 2(1 - uu) p ■ p (pj • p 2 + p 2 ■ pi) 

+ -P' P'{u 2 (l + u 2 )p 2 ■ p' 2 + u 2 (l + u 2 ) Pl -p'A. (A.6) 
uu L J 

Here and everywhere u = 1 — u and 

P = P!+p 2 , p' = pi + p' 2 • (A.7) 

The computation of the diagrams with one nonabelian vertex runs along the same line. E.g., 
for Fig. |3] (g) we have 

AT(aj,y) = -8^ £ duu 2 u J dp„p 2 _ J d 2 z d 2 z' tr [t a t b u(x)t c u ] (y)) if cde u ad (z)u be (z') 

f d 2 p l d 2 p[ d 2 p 2 d 2 p' 2 i( Pl+P2 ). a! -i(p' 1+ p^. y -i(p 1 -p' 1 ). z -i(p 2 - ? / 2 ). z ' 
7 (2vr) 2 (2tt) 2 (2vr) 2 (2vr) 2 

v D2(Pl,P2,P'l,P' 2 ) (AR s 

o/i ;\2 /- 9 o\ /- ;2 / 2\ ' 

Pl (Pl + P 2 ) OPl + «P 2 ) («Pl + WP2 J 

with 

D 2 (pi,p 2 , pi, Pa) = r w>(p'i> -Pi - P 2 )P2)^ (Pi + P 2 ) 4a (Pi) 4 M (Pi) (P 2 ) d i (Ps) • (A-9) 
It results into 

p 3 . (mm) 2 L> 2 = uup 1 - p 2 (pi 2 - p' 2 2 ) - 2wp x • pi p 2 • (pi + p 2 ) + 2up 2 - p' 2 p 1 ■ (pi + p' 2 ) . (A. 10) 

Other three diagrams of the same topology can be easily obtained by means of symmetry argu- 
ments. The final result is given in Eq. (|64|). 



B Fourier transformation. 

In this appendix we give technical details on the Fourier transformation. To this end we use two 
simple formulae, the Chisholm representation of the propagator and the (^-dimensional Euclidean 
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momentum integral, respectively, 

1 



daa m -V ap ', (B.i; 



p2m r(m) JO 
d d P Jp-z-Ap 2 



(2vr) rf (4vr) rf / 2 A rf / 2 ' 1 ' 

and the formula (0) as well. 

For the sake of defmiteness, consider the diagram |] (g) which exhibits all features. 

Mix, y) = -4 (a s N c ) 2 In — ^ du ( d 2 z d 2 z'M(x, z)M(z, z')Miz', y) (B.3) 

xb Jo J 

d 2 p l d 2 p[ d 2 p 2 d 2 p' 2 PiaP^faM^P^) 

(2tt) 2 ( 2 7t) 2 (2tt) 2 (27T) 2 ( pi + p' 2 ) 2 + Mp 2) + ^2 



with 



xe -,, Pl - p .,-«-, te ^,.«.fj_ e ,( Pl+R) .,-, W M»-v_J_ e .« P , +R »-„-. W+rt .. 

1 MP! MP? 



/a/? (pi , P 2 ) = uu 5 a/3 (pi 2 - p' 2 2 ) - 2m p' la (pi + p' 2 ) p + 2m p^ (pi + p 2 ) a . (B.4) 

The first letters of the Greek alphabet stand for 2D transverse space with metric 5 a p = —g^a — 
diag(l,l), a, (3, 7, ... = 1,2. The tensor structure factorizes and can be evaluated separately. 
Consider the first term in curly brackets. The integral over unprimed momenta gives, using Eqs. 

(OH) 

d 2 p l d 2 p 2 ip ^ (x _ z)+ip9 . (x _ zl) Pi a P 2 B , R ,x 

(2tt) 2 (2vr) 2e pUupl + upl) 1 J 

i \ 2 (X- Z) a (x - Z')p 



,2n J (x — z) 2 (u(x — z') 2 + u(x — z) 2 ) 

Note that this formulae is correct only for u < 1. However, this should not bother us since the 
boundary is not reached due to limits on the final state mass resulting into the cutoff (|75|). 
The primed momenta are integrated out with the formula 

d 2 p[ d 2 p' 2 ip i i . ly - z) -. ip i ) .( y -z>) fapjp'iM m 6) 

(2-) 2 (2vr) 2 (Pi+p' 2 ) 2 (npi 2 + np' 2 2 ) 1 ' ) 

- <&f* -€) 1: d i s: d i - p b h« - ^ + v ^ } ' 

where 

V(y) = u(z-y)+u(z' - y) . (B.7) 
The action of differential operators in f a g reduces to substitutions 
d 2 d 2 



dz 2 dz 



a 



-Ap |2u - 1 - ~ (v (2m - l)V(y) + 2uw (2 - 2')) • V(y)} 
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d 


( 9 


r) 
(J 


dz a 


\dzfj 


dz', 


d 


f d 


d 


dz', 


\dz a 


dz' 



-2up U a p - 2 P - (vV a (y) + u(z- z%) Vp{y)\ 



(B.8) 
Finally, we get 

d 2 p[ d 2 p' 2 ip ,. (y _ z) _ lp , y(y _ zl) faf3(p'l,P' 2 ) (v , q] 

(2^) 2 (2 7 r) 2 (pi + p> 2 f (up', 2 + up' 2 2 " 

i \ 2 1 



2nJ (z- z' f (u(z-y) 2 + u(z' - y) 2 ) 
x {uu{z - z') ■ (2V{y) + (1 - 2u){z - z'))^ - 2uVp(y) {z - z% - 2uV a {y) {z - z'),} . 

Note that 2V(y) + (1 — 2u) (z — z') = (z' — y) + (z — y) is u-independent. The second term in 
the curly brackets in (|B.3| ) is analyzed along the same line. 

Let us note that the calculation of the diagram |3] (a) does not present a difficulty either since 
Di from Eq. ( |A.6| ) is a 'square' of the Lorentz structures f a , ( |B.4j ). Namely, 

J) 2 _ (UU) D x = ■^rfap{Pl,P2)fap{p'l,P2) ■ (B.10) 



Thus, all the Fourier transforms factorize and are given by Eq. (|B.9|) . In the computation of the 
diagram |] (a) it is instructive to use the identity 

V(x) ■ V(y) = \ {V 2 (x) + V 2 (y) - (x - y) 2 } , (B.ll) 

in order to simplify the algebra. 
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